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Using recent results in the field of quantum chaos we derive explicit expressions for the time scale of 
decoherence induced by the system-environment entanglement. For a generic system-environment 
interaction and for a generic quantum chaotic system as environment, conditions are derived for 
energy eigenstates to be preferred states in the weak coupling regime. A simple model is introduced 
to numerically confirm our predictions. The results presented here may also help understanding the 
dynamics of quantum entanglement generation in chaotic quantum systems. 
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I. INTRODUCTION 

Real physical systems are never isolated from the sur- 
rounding world and, as a consequence, nonclassical cor- 
relations (entanglement) are established between the sys- 
tem and the environment. This process, which leads 
to decoherence, has a fundamental interest since it con- 
tributes to the understanding of the emergence of classi- 
cality in a world governed by the laws of quantum me- 
chanics 

Remarkably, different states may decohere at drasti- 
cally different rates, and a small fraction of them may be 
particularly stable under entangling interaction with the 
environment P, [2]. Such states are "preferred states" 
of the system under the influence of the environment. 
(They are also called "pointer states" , a name given for 
the states of pointers of measurement apparatus in the 
study of the measurement problem [l|.) This concept 
is important in understanding what states will naturally 
emerge from a quantum system subject to decoherence. 

As discussed by Paz and Zurek in an interesting paper 
[il, depending on the type and strength of the system- 
environment interaction, different preferred states may 
arise during the decoherence process. In particular they 
consider the case of an adiabatic environment modeled 
by a quantum scalar field and interacting weakly with a 
system through a given type of coupling. In such situa- 
tion they show that energy eigenstates are good preferred 
states and hence are the natural representation of the 
quantum system. Notice that the adiabatic environment 
does not change the population of energy eigenstates of 
the system, implying infinite relaxation time. 

In realistic situations the environment is often nona- 
diabatic, with finite relaxation time. In this case, the 
relation between the relaxation time and the decoher- 
ence time is crucial and, typically, only when the former 
is much longer than the latter, energy eigenstates can 
be good preferred states. While a rough estimate of the 



relaxation time can be obtained via Fermi's golden rule, 
for the decoherence time the situation is more complex. 
Indeed in case of generic type of coupling and generic 
environment, it is hardly possible to obtain estimates by 
employing the master-equation approach used in 

In this paper we propose an alternative approach which 
takes advantage of recent progress in the field of quantum 
chaos: namely random matrix theory and the so-called 
fidelity I, i i, 0, i which is a measure of the stability 
of the quantum motion under system perturbations. Us- 
ing this approach we can estimate - via the fidelity decay 
of the environment - the decoherence time of the system 
for weak, generic system-environment couplings and for a 
broad class of environments. In particular, we need not 
restrict ourselves to the Markovian regime (in contrast 
to master-equation approach) where the bath correlation 
decay is faster than decoherence. Moreover, as an im- 
portant result of ours, by modeling the environment by a 
chaotic quantum system, we determine a critical border 
in the coupling strength below which energy eigenstates 
are shown to be preferred states, while above this border, 
the relaxation time and decoherence time have the same 
scaling with the coupling strength and therefore no defi- 
nite statements can be made. We also present numerical 
results which confirm the above picture. 



II. GENERAL THEORY 

Let us consider a quantum system S with a discrete 
spectrum, weakly coupled to a second quantum system 
£ as its environment . The total Hamiltonian is: 



H = Hs + eHi 



H, 



(1) 



where Hs and Hs are the Hamiltonians of S and £ re- 
spectively, and eHj is the weak interaction Hamiltonian 
(e <C 1). The time evolution of the whole system is given 
by |^'5£(i)) e-*^*/'^ |*5£(0)). The initial state is set as 
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a product state |*s£(0)) = IV's(O)) |0f(O)). The reduced 
density matrix {a\p^'^ {t)\P) — {a[TT£p{t)\(3) is obtained 
by tracing over the environment. 

Consider first the case with initial state \il)s{^)) = 
I a), where |a) denotes an energy eigenstate of Hs with 
eigenenergy Ea- We define the projection operators 
\a){a\ ® Is and = T^p^c ® ^£ where Ig is 

the identity operator for the environment degrees of free- 
dom. The whole Hilbcrt space can then be decomposed 
into two orthogonal subspaces, leading to: 



)\€{t))+^\XT,{t)), (2) 



where e|xa(i)) = !^a\^ ssit)) ■ The small parameter e is 
introduced to account for the fact that, in case of weak 
coupling, the second term in Eq. ([2]) remains small inside 
some initial time interval (see below). The normalization 
of ssit)) in Eq. ^ is unity to the first order in e. 

A simple derivation shows that the evolution of the two 
terms in Eq. ([2]) is given by the coupled equations, 

'^J^\€{t)) = Ht^mt)) + e^e^^°*/^(a|i7,|xz.W>, (3) 



i^-^lVait)) = exp 



n 



where H'/J = eHio^+Hs, \ri^{t)) = exp {iH^t/h) |xa(t)), 
Hia = {a\Hi\a), and Ha = 3^aH^a- 

It is evident from Eq. ^ that e^(xalXa) gives the pop- 
ulation that has leaked to the subspace associated with 
In the case of weak coupling, this population leak- 
age is initially very small. More precisely we have that 
e^iXalXa) ^ 1 up to times t <^ te, where te is the re- 
laxation time of the system. Then the second term in 
Eq. ([3]) can be safely neglected and, as a result, the envi- 
ronment is in the state |0f (i)) ~ e~'*^£° /^|0£(O)) while 
the system remains in the eigenstate \a) with a definite 
phase evolution. 

Consider now as the initial state a superposition of 
energy eigenstates |V's(0)) = Y.a'^aW). As in Eq. ^ 
we can write: 



I* 



S£ 



C^\a)\ct>i{t))+e\x{i)), (5) 



where the term \x{t)) — Yl,a (^a\Xa{t)) contains now con- 
tributions of transitions between different energy eigen- 
states. Note that the first term on the right hand side 
of Eq. ([5]) may be highly entangled even when e\x{i)) is 
small. For t <^ te, the term e\x{t)) in Eq. (O can be ne- 
glected and, by tracing out the environment, the reduced 
density matrix of the system writes: 

p-^ = {a\TT£pitm c e-^^^--^^^'/^C^C*pfpo.{t), (6) 
where fi3a{t) = (0^(i)|'/>a(*)> satisfies 

with V = Hip-Hi^ = {(5\Hi\P)-{'^\Hi\a). (8) 
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FIG. 1: Time dependence of the elements of the reduced 
density matrix of the qubit system S. Here A'' = 200 and 
the coupling strength e = 5 x 10~* is much smaller than the 
perturbative border ep 0.04 numerically computed from 
Eq. ([21 . The upper thin solid curve gives the diagonal matrix 
element p22it) for an initial state with |7/)s(0)) = i2). The 
lower thick solid curve gives the off diagonal matrix element 
Pviif) for an initial state which is a product state in which 
the system S is in an energy superposition state while the 
environment is in a randomly chosen state. The dotted and 
dashed curves give the theoretical Gaussian decay Eq. (|10|) 
and the exponential decay Eq. (|12p . respectively. 



Significantly, the quantity //3a (t) is simply the "fidelity 
amplitude" of the environment associated with the two 
slightly different Hamiltonians H'^^J and {H^{/ + e-V). 
For non-constant V, this fidelity amplitude usually de- 
cays with time for a ^ /3, then p^'^ also decays and there- 
fore decoherence sets in. In case of constant V in Eq. ([S]), 
Eq. ([7|) gives \fpa(t)\ ~ 1, hence, there is no decoherence 
induced by the environment in the first order perturba- 
tion theory discussed above. Notice that Eqs. (6) and (7) 
become exact in the particular case in which the coupling 
Hi commutes with Hs- This case has been considered in 
% 

It turns out, from the above considerations, that for 
weak, but generic coupling, energy eigenstates of Hs play 
a special role in the sense that, only for energy superposi- 
tion states (not single eigenstates), the decoherence pro- 
cess is associated with the instability of the environment 
under perturbation (fidelity decay). 



III. A MODEL WITH THE ENVIRONMENT 
MODELLED BY A QUANTUM CHAOTIC 
SYSTEM 

Let us now turn to an explicit estimate of the decoher- 
ence time of the system. To this end we can directly apply 
recent results on fidelity decay [1,0,H,01 which, as shown 
below, allows us to estimate the decoherence time for a 
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generic type of system-environment interaction and for a 
broad class of environments. This contrasts the situation 
of the master-equation approach with which only par- 
ticular types of interaction and environment have been 
treated [3'| while extension to more general situations is 
mathematically difficult. 

Let us assume that the environment is modelled by 
a quantum chaotic system ^ld\ . (Analogous strategy can 
be applied when the environment has a regular or mixed- 
type phase space structure.) Fidelity decay in such sys- 
tems has been studied via semiclassical methods as well 
as with random matrix theory, both giving consistent 
results. Specifically, it turns out that for initial states 
chosen randomly, the fidelity has typically a Gaussian 
decay below a perturbative border Cp and an exponential 
decay above this border [1, H, 0]. If we model the en- 
vironment Hs by a single matrix of dimension N taken 
from the so-called Gaussian orthogonal ensemble (GOE) 
[ll|, then the border Ep can be explicitly estimated and 
is given by Q 

2^ep"KE^ - '^vA, (9) 

where V^^ is the average of |(n|F|n')p with n ^ n' and 
V given by Eq. ([5]). Here \n) denotes the eigenstates of 
H^^J , A is the mean level spacing of H^^f , and is the 
variance of (n|y|n). 

Below the perturbative border, e < e^, the fidelity am- 
plitude decays as [1| 

|/,„(^)|c.e-^-^^/2''^ (10) 

Then, the decoherence time Td, characterizing the decay 
of off-diagonal matrix elements [see Eq.®] is given by: 

Td V2h/{eay) (X e^^, e < Sp- (11) 

Notice that this dependence of Td on e coincides with 
the one derived in Ref. even though in our case we 
do not assume an adiabatic environment. For e > Cp, 
\ff3a{t)\ has an exponential decay Q: 

|//3o(t)| ^ e-"/2'' with r = 27re2v^/A, (12) 
and T<i~7^A/[7^e2T7^]cxe-^ e > ep. (13) 

We stress that, like Eq. (O, also Eqs. ([TT]) and US]) are 
valid only if the time scale under consideration is much 
less than te- 

The next issue is if, and under what conditions, te is 
sufficiently large so that significant decoherence may oc- 
cur for t <^ Te- If this is the case, then energy eigen- 
states are much more robust than their superposition 
states. Let j^g) be one eigenstate of Hg and {Hj^j) 
be the mean square of the non-diagonal matrix elements 
{a'\{iJ,£\Hj\fi£)\a) {a ^ a'). One may estimate te by 
using Fermi's golden rule, i.e., 

TE^l/Roc e-^, where R = 2ne^ p{H] ^^^) /h (14) 
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FIG. 2: Same as in Fig.[T]but for e = 0.01 which is still below 
but close to the perturbation border. It is seen that the nu- 
merically computed \pi2 \ begin to deviate from the predicted 
Gaussian decay. Moreover the population decay (thin solid 
curve) becomes appreciable. 



is the transition rate in Fermi's golden rule with the 
on-shell density-of-states p approximated by the average 
density of all possible final states for the whole system 

Therefore, below the perturbative border (e < ep) the 
decoherence time Td and the relaxation time te scale as 
e^^ and e^^, respectively. It follows that for small enough 
e, we have Td <^ te and hence energy eigenstates are 
good preferred states regardless of the form of the cou- 
pling. On the other hand, above the perturbative border, 
both Td in Eq. (fT3|) and te in Eq. (fT4|) scale with e~^. In 
particular, in cases of very small Cp, the two time scales 
can be comparable even at weak perturbation and 
hence energy eigenstates may not be preferred states . 

We will now introduce a simple dimensionless model 
which will also allow an explicit numerical evaluation of 
the different time scales. Let 5 be a qubit system with 
Hamiltonian Hs = Ea\a){a\, a — 1,2 and E2 — 
El = 1. The environment is modeled by a dimensional 
matrix in the GOE and the interaction Hj is taken in 
the form of a random matrix as well. Specifically, in a 
given arbitrary basis |fc), the matrix elements {k\H£\k') 
and {ak\Hi\a'k') are real random numbers distributed 
according to a Gaussian with unit variance. Moreover 
we set the Planck constant Ti = 1. 

We have numerically integrated the time-dependent 
Schrodinger equation for this model thus obtaining the 
matrix elements of the reduced density matrix p'"^. In 
Fig. [ll we show the numerical results for parameters 
N = 200 and e = 5 x 10""*. The perturbative border can 
be numerically computed from Eq. ^ and it is found to 
be ep ~ 0.04. It is clearly seen that for weak coupling the 
decay of Pi2 is well predicted by the Gaussian decay in 
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FIG. 3: Dependence of times Td and te on the coupling 
strength e. Full circles and empty squares represent the nu- 
merically computed and te, respectively. The solid line is 
given by the theoretical expression Eq. pif) and the dashed 
line is given by Eq. (|14|l . Notice the good agreement between 
numerical data and theoretical predictions. The arrow indi- 
cates the perturbative border ep which clearly separates the 
two different scaling behaviors of r^. 

Eq. pn]) . By contrast, the change in the diagonal matrix 
elements of the reduced density is negligible, as shown by 
the upper thin curve in Fig. [TJ 

Figured] is drawn for the same parameters of Fig.[l]but 
for a larger coupling strength, comparable to the pertur- 
bative border. One notices deviations from the Gaussian 
decay and also an appreciable population change (upper 
thin curve) . Interestingly, we found that deviations from 
the Gaussian decay always goes with an appreciable pop- 
ulation change. Indeed if, for example, we deliberately 
weaken those off-diagonal coupling terms (diagonal cou- 
pling terms not touched) that are responsible for the pop- 
ulation decay, then the Gaussian decay can be recovered 
while the population decay becomes negligible. 

Finally, if we further increase the coupling strength 
above the perturbative border, then the exponential de- 
cay of the off-diagonal matrix elements [Eq. ([T^] is in- 



deed observed (not shown here). 

We have carefully studied the scaling behavior of te 
and Td as functions of the coupling strength e and the 
results are shown in Fig. [31 Numerically is defined as 
the time scale over which decays by a factor of 1/e, 
and Te is defined as the reciprocal of the slope of P22 i^) 
in the initial interval of time in which Fermi's golden rule 
is valid. It is seen that numerical data nicely agree with 
analytical predictions. In particular, one can distinguish 
between the two scaling behaviors of t^ (e^^ and e^^) 
separated by the perturbative border ep. Below this bor- 
der, Td 'ti Te, implying that energy eigenstates are much 
more stable than energy superposition states. 

IV. DISCUSSIONS AND CONCLUSIONS 

Entanglement-induced decoherence within a closed to- 
tal system is also referred to as "intrinsic decoherence" 
[31 ■ The above analysis directly indicates the existence 
of preferred states of intrinsic decoherence. As such, the 
dynamics of entanglement generation between two sub- 
systems depends strongly on the coherence properties of 
the initial state. This is of interest to studies of entangle- 
ment generation in classically chaotic systems p^ . Fur- 
ther, our results might also shed light on the dynamics of 
quantum thermalization processes within a closed system 
where energy eigenstates also play a special role. 

In summary, for weak but generic coupling between 
two quantum subsystems, energy eigenstates are shown 
to play a special role in the entanglement-induced deco- 
herence process. The quality of the energy eigenstates 
as preferred states is also analyzed in terms of a simple 
dynamical model which allows for numerical analysis. 
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